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ABSTRACT 

We investigate an intensification mechanism for tlie magnetic field near the 
^ base of the solar convection zone that does not rely on differential rotation. Such 

pLf mechanism in addition to differential rotation has been suggested by studies of 

flux emergence, which typically require field strength in excess of those provided 
by differential rotation alone. We study here a process in which potential energy 
of the superadiabatically stratified convection zone is converted into magnetic 
energy. This mechanism, know as "explosion of magnetic flux tubes", has been 
previously studied in the thin flux tube approximation as well as two-dimensional 
MHD simulations, we expand the investigation to three-dimensional MHD sim- 
ulations. Our main result is that enough intensification can be achieved in a 
three-dimensional magnetic flux sheet as long as the spatial scale of the imposed 
perturbation normal to the magnetic field is sufficiently large. When this spa- 
tial scale is small, the flux sheet tends to rise toward the surface, resulting in a 
significant decrease of the magnetic field amplification. 

Subject headings: Sun: interior — Sun: dynamo — Stars: interiors 



Introduction 



The solar magnetic field is thought to be generated by dynamo action (jParkerl Il955l ). 

In most solar dynamo models differential rotatioi i is the key process that generates strong 

toroidal field near the base of the convection zone. ( Choudhuri et al]|l995 : Dikpati fc Charbonneau 
19991 : iHotta &: Yokoyamall2010al Jbl) . The mean energy density of dynamo-generated magnetic 
field is determined by the energy density of differential rotation as: 



— ~ —p(I\v] 
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where -B, p, and Af denote the strength of magnetic field, density and typical difference of 
the fiuid velocity at the tachocline. At the tachocline these values are p = 0.2 g cm~'^, Av = 
27r i?hnAr2 = 27r X 5 x 10^ ° cm x 50 nHz = 15700 cm s~^ (see the results of helioseismology 
in [Thompson et al.l l2003l ). where i?bc and AQ are the radius at the tachocline and the 



difference of angular velocity between the top and the bottom of the tachocline, respectively. 
In addition, if we assume a linear velocity profile and formally integrate that over the width 
of the shear layer, we have a drop by another factor of 3 in energy. Using these values the 
dynamo-generated magnetic field is estimated to be approximately i? ~ 1.4 x 10'^ G. This 
estimate is similar to the values found in mean-field models with Lorentz force feedback 



by iRempell (120061 ) . Note that in iRempell (120061 ) latitudinal shear mostly generates toroidal 
magnetic field and that differential rotation is continuously replenished through the A-effect 
(parameterized turbulent angular momentum transport), which explains why 10^ G can be 
reached with only small feedback on differential rotation. 

On the other hand, theoretical studies of flux emergence and their comparison with 
observed sunspot properties (Joy's law, low latitude emergence ) imply substantially stronger 
fleld near the base of the convection zone. I Weber et al.l (120111 ) suggest that a toroidal fleld 
strength of about 50 kG is required at the base of the convection zone in order to reproduce 
th e statistica l properties of tilted sunspot pairs in rising flux tube simulations (see also review 
byHliooJ. 



Above results suggest that additional mechanisms for the ampliflcation other than ro- 
tational shear are req uired. One ca i ididate of such a mechanism is the "explosion" of risiri g 
magnetic flux tubes. (jParkeiill994j : iMoreno-Insertis et al.lll995l : iRempel fc Schusslerll200ll ). 
It is summarized as follows: Suppose that there is a superadiabatically stratifled atmosphere 
like the solar convection zone and a dynamo-generated toroidal magnetic fleld. When the 
gravity has the proflle of g{r) oc r~^, where r denotes the radius, the distribution of gas 
pressure Pe{r) is expressed as 



Pc{r) = Pba 



1 + V 



base 
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using a constant temperature gradient as 



V = 



dlogT 
dlogp 
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Pbase and ifbase dcuotc the gas pressure and the pressure scale height at r = rbase, respectively, 
where rbase is the location of the base of the convection zone. As the toroidal fleld buoyantly 
rises as a flux tube, the thermal property in the tube vary adiabatically while the surrounding 
gas in the convection zone is stratifled in a superadiabatic manner. The thermal exchange 
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is mainly by radiative diffusion in the solar convection zone and the time scale of it is at 
least several tens of years for a 1000 km flux sheet, which is much longer than the time scale 
of flux emergence on the order of months. Therefore the internal distribution of the gas 
pressure is expressed as: 



1 + V 



ad 



base 



- 1 



1 1/Vad 



(4) 



where Vad denotes the adiabatic temperature gradient. We note that V > Vad and the 
pressure of the internal medium of magnetic flux decreases slower with r than that of the 
external medium. Suppose that the pressure is in balance at the base of convection zone as: 



Pil'^ba 



+ 



8tt 
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As the decrease of the pressure in the tube is slower than that of the surrounding plasma, 
these pressures become the same at a certain height r = rexpi, i.e., Pi(rexpi) = Pe('^expi)- This 
height is called the explosion height. Even after passing the explosion height, the flux tube 
continues to rise up to the surface of the sun r = rt. During this process, mass in the flux 
tube is swept out due to the imbalance betw 'een the pressure gradient and the gravitational 
force along the magnetic field (see Fig. 2 of iRempel fc Schusslerll200ll ). We specify the gas 
pressure distribution Pi^t^r), after this sweeping of mass. Since the gas pressure is balanced 
at the surface, i.e. Pit('"t) =Pe{rt), the pressure distribution is expressed as: 



Pi,t{^) = Pcin) 
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where a subscript "t" denotes the value at the surface. If both ends of the the rising flux tube 
stay at the initial position (r = rtase), the amplified magnetic field has the energy density of 



td2 

amp 

Stt 



Pc('"basc) ~ Pi, t ('"base); 



(7) 



Since the pressure is in balance at the surface with the amplified field, the magnetic field 
is amplified up to the strength for which the explosion height approaches the surface. The 
asymptotic field strength is independent of the initial magnetic field and given by the su- 
peradiabaticity of the convection zone. 



Moreno- Insertis et al.l (119951) iny e stigat e explosion using one- dimensional thin flux tube 



approximation. iRempel fc Schiisslen (120011 ) show using two-dimensional MHD calculation 
that magnetic field can be amplified up to 10^ G in the solar convection zone. The time 
scale of amplification is determined by the initial magnetic field and is < 0.5 year when the 
initial plasma beta is /3 < 10^. 
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In this study, we investigate the possibihty of an explosion of a magnetic flux and its 
consequential intensification by MHD simulations in three-dimensions, i.e. in more realistic 
situations, which have not yet been investigated in previous studies. 



Model 



The three-dimensional magnetohydrodynamic equations are solved in Cartesian coor- 
dinates {x,y,z), where x and y denote the horizontal directions and z denotes the vertical 
direction. Equations are expressed as: 
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where p, p, v, B, g and e denote density, pressure, fluid velocity, magnetic fleld, gravitational 
acceleration, and total energy per unit volume. The value of a speciflc heat ration is set as 
7 = 5/3. Note that since we use the equation of state for the perfect gas, the internal energy 
per unit volume is expressed as p/(7 — 1) in eq. ffT^ 



m 



The divergence fre e condition, i.e. V-B = 0, is maintained using the method introduced 
Dedner et al.l ( 12002| ). The simulation domain is (—7.5, —20, —1) < {x j H\,,y j H\„ z / H\,) < 



(7.5, 20, 3) and the total grid number is 256 x 128 x 384. We use a forth-order centered 
finite difference scheme for spatial-derivatives in combination with a forth-order Runge- 
Kutta integration in time. We do not use any explicit viscosit y or magnetic diffu sivity, 
but apply artificial diffusivities to maintain numerical stability (IRempel et al.ll2009l ). We 
use periodic boundary conditions for all variables in horizontal directions. At the top and 
bottom boundary, stress-free and non-penetrating boundary conditions are adopted, i.e. 
dvx/dz = dvy/dz = and = 0. The values of density and pressure at these boundaries 
are derived from solving eq. ( 120|) -( 123|) . At the top boundary, the magnetic field is set as 
Bx = By = 0, and dBz/dz = 0. At the bottom boundary the magnetic field is set as, 
dBjdz = dBy/dz = dBjdz = 0. 
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The profile of gravitational acceleration is assumed to be 

2 



9h 



rh 
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where rb = 5i^b is the location of 2; = from the center of the gravity source. We note that 
the subscript 'b' means the value at 2; = 0. Here Hb is used instead of i^base since the setup 
of the simulations here is similar to but not identical with the model in Section [T] The initial 
magnetic field is set as a flux sheet with a perturbation as: 

{z-z,r 



Bx = Bof exp 
a = ao + ai exp 



/ = 
Bo 



Zo + zi exp 

Oo 

a 



X 

dx 

X 



y_ 

dy 

y_ 



Svrpb 



(14) 
(15) 

(16) 
(17) 

(18) 
(19) 



where a, and Zc denote the thickness and the center of a flux sheet, respectively, dx and 
dy are the range of perturbation in x and y direction, respectively. We specify ao = 0.2ifb, 
Oi = 0.25ifb5 zq = O.li^b, Zi = 0.25ifb; /3 = 300 and dx = 2ifb- dy is a free parameter 
from 2ifb to SH^,. As shown above, we are considering a toroidal flux sheet as a source 
of the magnetic flux prior to the intensification because the magnetic field stretched by 
the f2-effect at the bottorn of the convection zone is suggested to have a sheet structure 
according to iRempell (120061 ). Then is calculated to satisfy the divergence free condition. 



i.e. dBx/dx + dBz/dz = with i?^ = at the all boundaries and By = in all the domain. 
In order to mimic the situ ation of the solar convection zone, we adopt a similar approach 



as lRempel fc Schiissleii (I2OOII ). The upper layer {z > 0) is adiabatically stratified (convection 
zone: 6 = 0) and the lower layer {z < 0) is subadiabatically stratified (radiative zone: 
6 = —0.2), where S is the superadiabaticity. The magnetic layer is placed at the interface 
between both regions. Since our convection zone is adiabatically stratified we add additional 
entropy within the flux sheet to allow for a buoyant rise and subsequent explosion. Using 
this setting, the pressure scale height in the flux sheet is longer than that of the external 
medium and the explosion occurs around 2;expi ~ -f^bCp//3As, where As is the added entropy. 
Using an adiabatic stratification in the upper layer allows us to study the explosion process 
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in a relatively simple situation without thermal convection; this leads, however, also to 
limitations since some key processes are not catured in our setup. Beside not capturing 
the influence from convection, our magnetic flux sheet is not stored in stable layer having 
neutral magnetic buoyancy, i.e. overshoot region at the base of the solar convection zone. 
The flux sheet is however stabilized against buoyancy to some degree by lying on top of our 
stably stratified lower layer. We will investigate the explosion process in a superadiabatic 
convection zone with fully developed thermal convection in a future study. 



exp 



[Z - Zr 



The initial conditions for pressure and density are derived from solving the equations: 

''^ (20) 
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We specify As = O.OScv. Using this parameter the explosion height is around Zexpi = O.lSifb- 
We do not yet aim at a realistic simulation with solar parameters but rather consider a model 
that allows us to study the explosion and the following intensification in (artificial) isolation. 
If there is a superadibatically stratified convection zone, turbulent thermal convection is 
generated and the situation becomes more complicated. 



3. Result 

In this study, we take the y-component of the perturbation on the initial magnetic field 
{dy) as a free parameter. Fig. [1] shows the strength and the field lines of the magnetic field for 
simulations with dy = 2H\^ and dy = QH\^. We present here only the solutions in the center 
region cut at ?/ = to enhance the clarity of the presentation. The magnetic layer fills in the 
horizontal directions the whole domain outlined by the white lines. In our three-dimensional 
situation a value of dy = 6ifb (Fig. [T]i) leads to a substantial intens ification of magnetic 



field, which is similar to the two-dimensional calculation reported by iRempel fc Schiissler 



(1200 ll ). At the foot-points the magnetic field lines remain almost parallel to the x-axis for 
dy = QH^. On the other hand, the amplification of magnetic field is significantly less with 
dy = 2H\y. In the result with dy = 2ifb, the foot-point magnetic field rises upward and the 
magnetic field lines at the back {y < 0) are bent to the center {y ~ 0) of the calculation 
domain. At this stage, it is suggested that the intensification of magnetic field by explosion 
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process is more likely to occur with a longer perturbation, which is perpendicular to the 
main magnetic field. 

Fig. [2^ shows the temporal evolution of magnetic energy averaged over the y = plane. 
Around t = 120i7b/cb, where Cb is the speed of sound at 2; = 0, the initial intensification of 
magnetic field with explosion process is completed in all cases, i.e. the top of the rising fiux 
is reaching the upper boundary. Fig. [2b shows the dependence on dy of the mean magnetic 
energy at t = 120ifb/cb- It is qualitatively clear that a larger scale dy causes a more effective 
intensification. 

In the following discussion, we compare the results with different dy at the time when 
the apex of the rising magnetic flux sheet reaches z = 2ifb in each case. Fig. |2h shows 
the proflle of the y-component of the fluid velocity averaged over — O.li^b < a; < O.lifb and 
< z < 2H\, at ?/ = 0. Fig[21i shows the dependence of the maximum value of the averaged 
y-component fluid velocity. The flow parallel to the |/-axis is generated by a pressure gradient 
dp/dy. As the magnetic flux rises, the gas pressure in the central region decreases and it 
causes a converging flow in the y-direction. The amplitude of the fluid velocity is estimated 
as Vy = rAp/ (pdy), where Ap denotes the perturbations of the gas pressure with respect to 
the ambient medium and r the time intervall over which the acceleration takes place. For 
our estimate we use values of p ~ pb, and r ~ lOOifb/cb, The difference of the pressure is 
estimated to be Ap ~ O.OOlpb which is the same order of Ph/ f3, where pb and pb are the 
density and pressure at z = 0. Thus we obtain Vy = 0.1H\^C]^/dy. The dashed line in Fig. 
|2Ji corresponds to this equation. Overall the result with dy = 4:Hy,, 6H]j, and 8ifb are well 
reproduced. The result with dy = 2iJb is slightly smaller than the value estimated from 
Vy = 0.1/ dy. The reason is that the center region is already fllled with mass, resulting in a 
smaller Ap. In addition, a small value of r with dy = 2H\^ causes the deviation. Fig. |3^ 
shows the value of f perp along the fleld lines, where f perp is the fluid velocity perpendicular to 
the magnetic fleld in the x-z plane. A positive (negative) value of f perp indicates a flow which 
moves the magnetic flux upward (downward). Using small (large) dy, the foot-point magnetic 
flux move upward (downward). Fig[3]D shows the velocity parallel to the magnetic fleld along 
magnetic fleld lines. The larger velocity is generated with larger dy. These correspond to 
the ampliflcation of the magnetic fleld. 



4. Discussion 



Our results show that the efficiency of the magnetic intensification is dependent on 
the spatial scale of the perturbations {d,/]. The r eason is understood as follows: In a two- 



dimensional situation (IRempel &: Schiisslerl I2OOII ) or using large value of dy{> Gifb), the 
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pressure gradient dp/dx leads to strong converging flows toward the center of the rising flux 
sheet within the y = 0-plane. The flow component parallel to the field leads to amplification, 
the flow component perpendicular to the field prevents the buoyant motion around the legs of 
the rising flux and consequently plays a role for pinning down the foot-points (Fig. H^). On 
the other hand, using small dy{< 2ifb), a converging motion in the ^/-direction is promoted 
by the large pressure gradient dp/dy and the center of the explosion region is filled with the 
mass. Thus the pressure around the foot-point becomes large and the flow is suppressed. 
Consequently magnetic flux moves upward and is less amplified (Fig. Wp). If the foot- 
point magnetic field is fixed at the initial position, the amplification of magnetic energy is 
estimated with eq. ([7]). In a three-dimensional situation, however there is a possibility that 
the foot-point magnetic flux moves upward, thus the amplified magnetic energy is estimated 
using eq. ([7]) at r = rf instead of r = rbase, where rf denotes the location of the foot- 
point magnetic flux. Fig. O shows the dependence of the amplified magnetic energy by the 
explosion process on the location of foot-point with the parameters of our calculation. If the 
foot-point magnetic field is fixed in the initial position, the magnetic energy of ~ 2 x 10~^pb 
can be generated by the explosion process. The result with large dy{= GH^,) well reproduces 
this estimation (see Fig. [T]i). Note that Fig. ^ shows the mean magnetic energy in the 
y = plane, leading to a lower values of ~ 6 x 10~^ph at maximum. Fig. [5] shows that the 
efficiency of the intensification significantly decrease as the foot-point magnetic field rises. 
It is confirmed that even using larger /3(= 600), we can obtain similar results, i.e. strong 
amplification is achieved with larger dy{> 6Hi,). 



5. Summary 

We investigate an intensification mechanism of magnetic field near the base of solar 
convection zone following the "explosion" process using three-dimensional MHD simulation. 
Such a mechanism in addition to differential rotation has been suggested by studies of flux 
emergence, which typically require field strength in excess of those provided by differential 
rotation alone. Our main result is that the enough intensification can be achieved even in a 
three-dimensional situation as long as the spatial scale of the imposed perturbation in the 
direction normal to the magnetic flux sheet is sufficiently large. When this spatial scale is 
small, the efficiency of intensification of magnetic field significantly decreases. Here a strong 
converging flow in the direction perpendicular to the field (y-direction) is driven, which leads 
to weaker flows in the x-direction. The latter are crucial for the field amplifiction as well as 
stabilization of the foot-points in our setup. 



In this study, we used a similar setup as iRempel Sz Schiissler 



(I2OOII ). which allows to 
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study the explosion process in an adiabatically stratified atmosphere, avoiding the complex 
situation caused by presence of turbulent thermal convection. Using a more realistic setup 
with a superadiabatic convection zone and storage of magnetic flux in a subadiabatic region 
beneath could influence our results in the following way: 

If the initial magnetic flux is stored in a subadiabatic stratification the foot-point mag- 
netic field is more stable and the criterion for the length scale of the perturbation might 
be relaxed, i.e. also shorter perturbations could lead to stronger amplification. However, 
even with more realistic storage in a subadiabatic region, there is no force that could stop 
the horizontal motion of the foot-points (buoyancy only works in vertical direction). As a 
consequence also here a fundamental difference between tube-like and sheet-like behavior 
will persist. 

The presence of turbulent thermal convection could significantly influence our results 
since during the explosion process the strength of magnetic field near the top of the rising 
flux drops significantly. Thus, it is possible that the connection between top and foot-point 
magnetic field is disturbed by turbulence, resulting in a reduction of field intensification. 

Therefore more realistic simulations including a stably stratified overshoot region/radiation 
zone and a superadiabatic convection zone are required in order to investigate the detailed 
conditions under which the explosion process leads to substantial intensification of magnetic 
field. 
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(b) d =6Hu: t=25HJc 




Fig. 1. — The volume-rendered magnetic energy and the magnetic field lines. The panels a 
and c (b and d) show the results with dy = 2H\^ {dy = 6b). For convenience of presentation, 
we show only the solution in the center region cut at y = 0. The color shows the magnetic 
energy normalized b y the gas pressure at z = (ph). This y isuahzation is created using the 
software of VAPOR jcivne fc Ra^^boosi Icivne et al.l 120071 ). 
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Fig. 2. — (a) Temporal evolution of the mean magnetic energy |i?^/87r| on the y = plane, 
(b) The dependence of the mean magnetic energy at t = 120iJb/cb on dy. The energy of 
magnetic field is normalized by lO^^pb- (c) The profile of y-component fluid velocity averaged 
over — O.lifb < X < O.lifb and < z < 2H\y at ?/ = 0. (d) The dependence of the maximum 
value of the averaged y-component velocity on dy. The dashed line denotes the dependence 

of Vy = 0.1 /dy. 
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Fig. 3. — (a) Averaged fperp along magnetic field lines, where Vperp is the fluid speed perpen- 
dicular to the magnetic field. We choose the magnetic field lines which start x = y = and 
z = 0, 0.25ifb- (b) Averaged velocity parallel to magnetic field fpara is shown. The same 
magnetic field lines for the panel (a) are choosen. 




Fig. 4. — Schematic picture of the explosion process in three-dimensional situation with (a) 
large and (b) small dy. 
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Fig. 5. — The dependence of the amphfied magnetic energy by the explosion process on the 
location of foot-point. 



